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ABSTRACT 

We  consider,  in  the  effective  mass  approximation,  the  problem 
of  a  hydrogenic  atom  in  a  strong  static  magnetic  field.  The  ground 
state  and  the  lowest  lying  states  of  angular  momentum  ±fi  are  cal¬ 
culated  variationally,  first  fora  simple  parabolic  conduction  band 
and  then  for  the  nonparabolic  conduction  band  of  InSb.  Both  the 
use  of  improved  variational  trial  functions  and  the  incorporation 
into  the  theory  of  the  nonparabolic  conduction  band  are  novel  fea¬ 
tures  of  the  present  calculation.  The  results  are  compared  with 
previous  work.  Nonparabolic  effects  are  found  to  be  significant 
for  the  energy  of  the  lowest  lying  state  of  angular  momentum  +  ft. 

The  weak-coupling  polaron  theory  is  used  to  derive  an  upper  bound 
to  the  shift  of  the  ionization  energy  of  the  ground  state  of  the  im¬ 
purity  due  to  electron- LO  phonon  interaction.  Polaron  effects  on 
the  ground  state  ionization  energy  are  thereby  shown  to  be  negli¬ 
gible  in  InSb.  The  bound  obtained  may,  however,  prove  useful 
in  estimating  polaron  effects  in  materials  with  larger  coupling 
constants. 
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SHALLOW  DONOR  LEVELS  OF  InSb  IN  A  MAGNETIC  FIELD 


I.  INTRODUCTION 

1  2 

Cyclotron  resonance  and  Hall  coefficient  measurements  In  “pure"  n-type  InSb  at  low 
temperatures  indicate  that  free  carriers  become  bound  to  ionized  donors  in  moderate  to  strong 
magnetic  fields. 

This  so-called  magnetic  freeze-out  of  conduction  band  electrons  had  been  predicted  by 
3 

Yaf^t,  Keyes  and  Adams  (YKA),  who  estimated  that  the  application  of  a  magnetic  field  of  25  kG 
lO  a  hydrogenic  donor  in  InSb  would  increase  its  binding  energy  by  approximately  a  factor  of  four. 

4 

Wallis  and  Bowlden  (WB)  extended  the  ground  state  calculation  of  YKA  for  a  hydrogen  atom 
in  a  strong  magnetic  field  to  a  variety  of  excited  states  and  showed  that  strong  optical  transitions 
should  occur  between  the  ground  state  and  the  lowest  lying  states  of  angular  momentum  Suer 

transitions  have  recently  been  observed  and  their  energies  measured  as  a  function  of  magnetic 
field5 

We  can  now  inquire  to  what  extent  the  model  of  isolated  hydrogenic  impurities  treated  in  the 
effective  mass  approximation  can  account  quantitatively  for  the  optical  and  Hall  effect  data  in 
InSb.  The  purpose  of  the  present  work  is  to  improve  previous  calculations  of  donor  energy  levels 
by  introducing  the  following  refinements:  (1)  improved  variational  trial  functions,  (2)  accounting 
for  the  nonparabolic  nature  of  the  conduction  band,  (3)  estimation  of  the  effect  of  electron- LO 
phonon  interaction. 


II.  VARIATIONAL  TREATMENT  FOR  PARABOLIC  BAND 

The  energy  levels  of  a  hydrogenic  impurity  in  a  simple  parabolic  band  subjected  to  a  uniform 
magnetic  field,  H,  are  eigenvalues  of  the  Hamiltonian  (1): 

ft  =  ^  +  ~Y  (xpy  -  ypx>  +  I  muc2(x2  +  y2)  +  XV(r)  ,  (1) 

where 


V(r) 


c  ,  2  x  2  2,1/2 

£q(x  +  y  +  z  ) 


2 

e _ 

O 


(2) 


is  the  static  dielectric  constant,  m  is  the  band  mass,  cjc  is  the  cyclotron  frequency,  eH/mc, 
and  we  have  used  the  vector  potential  A  defined  by  A  =  H/2(— y,  x,  0),  corresponding  to  a  uniform 
magnetic  field  in  the  z-direction.  We  take  the  parameter  X  equal  to  1. 
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1  _  I 

If  our  unit  of  length  is  (ft/2mto^)2  =  2  2  x  cyclotron  radius,  and  our  unit  of  energy,  we 

can  rewrite  (1)  in  dimensionless  form  in  cylindrical  polar  coordinates  as 


ft  =  H/ficj 


*0  +  Ki 


=  -  (if  +  i  1  t  1 

°  \dpz  »  *  p2 


d(p 


-(,/2  r‘/2/(p2*z2)1/2 


3z 


(3) 


where  L  is  the  angular  momentum  operator  l/i  0/ £)</?  with  integer  eigenvalues  M, 


y  = 


tied 
_ c 

2R 


and  R  is  the  hydrogenic  Rydberg. 
The  Schrodinger  equation 

K'd  =  Eil 


(4) 


is  not  completely  separable  and  we  are  forced  to  resort  to  a  variational  solution.  YKA  have 
pointed  out  that  in  the  limit  y  °°  (strong  field  limit)  the  eigenfunctions  of  H  take  the  product 
form 


f(p)  g(z) 


(5) 


1 


and  in  fact  f(p)  =  <pNM(p)  where  KQ(/?NM  =  (N  +  y) 

Since  <p  =  e  p  '  8  and  <P  ^  ^  =  p  e  ^  ^  e^,  one  is  motivated  to  use  for  the  ground  state 

V-i' 

trial  function 


2/q  2 
-P  /8a 


-z2/8b2 


(6a) 


and  for  the  lowest  lying  M  =  ±1  states: 


P  e 


2/fl  2  2  /qk 2 
-p  /8a  -z  /8b 


1<P 


-p2/  8a2  -z2/8b2  -i<p 

pc'  e  e 


(6b) 


(6c) 


2  2 

YKA  actually  took  (6a)  as  their  ground  state  trial  function,  varying  both  a  and  b  ,  while  \VB 

2  2 

used  (6a-c)  as  trial  functions  fixing  a  =1  and  varying  b  to  obtain  lowest  energy.  We  remark 
that  since  M  is  a  good  quantum  number  of  ft,  a  trial  function  which  is  an  eigenfunction  of  I , 
with  eigenvalue  Mr  yields  an  energy  which  is  an  upper  bound  to  the  energy  of  the  lowest-lying 
M  =  M'  eigenstate  of  ft. 

The  wave  functions  (6a-c)  have  two  important  weaknesses  as  variational  trial  functions: 

(1)  the  H  =  0  eigenfunctions  of  ft,  being  hydrogenic  wave  functions,  can  not  be  written  in  the  prod¬ 
uct  form  (5)  for  any  choice  of  f  and  g;  thus  any  trial  function  of  form  (5)  is  likely  to  get  pro¬ 
gressively  worse  as  H  decreases;  (2)  even  at  high  fields,  wave  functions  (6a-c)  decay  too  rapidly 
with  increasing  z  at  constant  p  (because  the  effective  potential  for  motion  in  the  z-directionb 
vanishes  as  z  -*  °°). 
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To  eliminate  these  shortcomings  we  have'  chosen  the  following  trial  functions: 


P  e 


2  /  2  .  2  2.1/2 

-p  /8a  c~k(p  +OZ  ) 

2/o  2  /  2  ,  2  .  _  2. 1  /2 

-p  /8a  -<(p  +  Qfz  +/?  )  ±i<p 


M  =  0 

M  =  ±1 


(7a) 

(7h) 


where  a  .  k.  rv  and  (3  are  variational  parameters  which  are  adjusted  to  optimize  independently 
each  of  the  above  functions.  Note  that  (7a)  and  (7b)  become  Is  and  2p  hydrogenic  eigenfunctions 
respectively  as  a  -*  oo>  a  -*  1  and  /3  —  0.  Further,  they  decay  in  a  realistic  way  as  z  -**  00  with 
p  fixed.  The  price  we  pay  for  our  improved  functions  is  that  the  required  integrations  can  no 
longer  be  performed  analytically. 

In  Table  1  we  tabulate  ionization  energies,*  deduced  from  the  variational  calculations  of 
WB,  YKA  (wave  functions  (6a-c)  with  a  and  b  varied)  and  the  present  calculation  denoted  PC. 
All  energies  are  tabulated  in  units  of  the  hydrogenic  Rydberg. 


TABLE 

1 

y 

Ei>M=0) 

eyka(m=0) 

E^)I(l(M=0) 

Ewk(M=1) 

4ka(m=1) 

Ep^UVI  =  1 ) 

0 

0.000 

0.8488 

1.000 

0.000 

0.2264 

0.2500 

i 

1 .346 

1,523 

1.661 

0.8564 

0.8919 

0.9118 

5 

2.520 

2.615 

2.750 

1.664 

1.686 

1.717 

25 

4.515 

4.565 

4.733 

3.100 

3.113 

3.183 

100 

7.  174 

7.200 

7.457 

5.095 

5.1030 

5.261 

Table  1  indicates  that  the  wave  functions  (7)  give  lower  energies  than  (6)  over  the  range 
of  y  calculated.  Notice,  however,  that  the  percentage  deviation  of  EyKA  ^rom  ^pc  decreases 
with  increasing  field  and  is  smaller  at  given  field  for  the  M  =  1  level  than  for  IV1  =  0.  This  last 
result  might  have  been  anticipated  from  the  fact  that  the  M  =  1  wave  function  vanishes  at  the 
Coulomb  center.  This  means  that  the  effect  of  the  Coulomb  field  is  weaker  on  the  M  =  1  state 
than  on  the  M  =  0.  The  strong  field  product  wave  function  (5)  should  therefore  be,  at  the  same 
value  of  H,  a  better  approximation  to  the  M  =  1  state  than  to  the  M  =  0  state. 

Because  it  can  be  deduced  directly  from  cyclotron  resonance  experiments,  the  quantity 
A(H)  =  E^(M  =  0)  —  E^(M  =  1)  is  of  particular  interest.  From  Table  1,  we  calculate  that  A® 

/  T  T  \  * 

is  always  greater  than  and  the  percentage  difference  becomes  particularly  marked  at 

lower  fields. 

Finally,  we  expect,  as  noted  above,  that  our  trial  functions  (7a,  b)  should  also  be  quite 
satisfactory  at  y  £  1,  where  product  functions  of  the  form  (5)  are  suspect.  Since  there  is  some 


*  The  ionization  energy  of  an  impurity  level  is  defined  as  follows.  Let  E^  ^  be  the  energy  of 

’  o 

the  ith  eigenstate  of  (1)  with  X  =  XQ.  Then  the  ionization  energy  of  this  eigenstate  is 

lim  (E.  —  E.  .).  Thus,  in  the  parabolic  band  case  for  M  =  0  and  M  =  —  1  states  the  ionization 

i,  x  1,  i 

A  (I)  0  0 

energy  in  Rydbergs  is  defined  by  E  (M  -  )  =  y  —  E(M  =  .  )  while  for  the  M  =  1  level  the 

m  1  1 

ionization  energy  is  E  (M  =  1)  =  3y  —  E(M  =  1). 
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experimental  interest  in  this  region  of  y,  we  tabulate  below  our  variational  results  for  the  hy- 
drogenic  ground  state  ionization  energy  with  y  <  3. 


TABLE  2 


y 

0.3 

0.5 

0.7 

1.0 

1.5 

2.0 

2.5 

3.0 

F(I) 

PC 

1.26 

1.39 

1.51 

1 . 66 

1.87 

2.04 

2.19 

2.32 

HI.  NONPARABOLIC  BAND 

In  order  to  study  the  effects  on  binding  energy  of  the  nonparabolic  structure  of  the  InSb  con¬ 
duction  bands,  we  turn  to  the  model  of  Bowers  and  Yafet.  This  model  gives  a  simple  and  accu¬ 
rate  description  of  the  conduction  band  Landau  level  energies  observed  in  InSb. 

BY 

The  BY  model  is  defined  by  the  8X8  Hamiltonian  matrix  H  f  given  in  the  position  repre¬ 
sentation  by 


H 


BY 


v  + 


<p+  lA)2 


2m 


<5  .  +  (p  +  -  A)  *  p  . 

nn*  p  c 


(8) 


where  fn,  is  the  energy  of  the  edge  of  band  nf  in  the  absence  of  magnetic  fields,  m^  is  the  elec¬ 
tron  mass  in  vacuum,  and  p.  ,  is  the  matrix  element  of  l/m  fi/i  v  in  the  band  edge  functions 

jnn*  v  j  e8 

u^  and  u^,.  We  use  the  same  band-edge  functions  and  notation  as  Pidgeon  and  Brown,  define 


our  impurity  Hamiltonian  by  the  matrix 

—  BY 

H  ,  =  H  ,  +  V(r)  6  , 

nn'  nn’  nn' 

and  look  for  conduction  band  solutions  to  the  eigenvalue  problem 


(9) 


y  h  ,f  =  Ef 
^  nn*  n’  n 

n'  =  1 


(10) 


The  functions  f  in  (10)  are  the  usual  envelope  functions  associated  with  the  total  electronic 


wave  function  given  by 
*  =  V 


f  u 
n  n 


Taking  our  zero  of  energy  at  the  valence  band  edge  in  the  absence  of  magnetic  field  we 
write  out  (10)  explicitly: 

k 

VT 

k 


<eg  -X')  +  P 

(Eg  -  A’)  f2  +  p 


K  f  K*  .  ,  kz  f  .  /  2  ,  f  x  K*.  \ 

3  V6  5  «^3  8  Z  6  V3  7/ 

K*,  ,  K  t  .  .  .  ( 2  .  ,  K  r  \ 

fA  +  -  f7  1  l - k  +  -  fR 

6  »ji  7  z  5  8/ 


K 

\T6 


=  0 


=  0 


Vr3‘^f. 

\'f.  =  p— f? 

4  \Tz  4 


( 1 1  a) 


(lib) 


(11c) 

(lid) 
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x.f  ap(JLf  )  . 

5  V-/6  1  46  Z> 

X'f,  =  p(—  k  f,  +  —  f,) 

6  2  1  */6  2' 

(A  +  x')  f ,  =  p(—  Kf,  +  -2-  f  ) 
7  Vx/3  1  nH  2/ 


<A  +  *  ’  f§ "  Pv^  1  41 


41 

/k  .  v 

(—  f,  +  —  K*f  ) 

\J1  1  ji  2/ 


(He) 

(Ilf) 

dig) 

(llh) 


where  X'  =  E  —  V(r),  K  =  (px  +  e/c  A^)  +  Kp^.  +  e/c  A^),  =  pz>  is  the  band  gap,  A  is  the 

spin-orbit  splitting  of  the  valence  band  and  for  all  values  of  j,  n  and  n*  for  which  p .  does  not 

vanish,  p  ,  =  p. 

*  *3  nn'  p 

Left  multiplying  (11a)  and  (lib)  by  X'  and  inserting  (llc-h)  into  (11a)  and  (lib)  gives  a  pair 
of  coupled  equations  for  f^  and  f^: 


X'(Er-X')  ^2 

- £ -  + 


P_  Ik 

Er,  \2 


KK*  +  |  K*K  +  f  k  2  +  X'K 


3  z 


11  1 
7  a  r-n  K  +  4  x’k  -  —  k 
3  A  +  X'  3  z  A  +  X’  z 


+  I  [*’.  K]  K*  +  i  124  K*}  K  +  - 

2  6  X1  3 


fX’.kl 

— T7~  k 
X’  z 


V  f  ^2 1  Ax'-K*i  i- 

j  *1  P  3  i - X1 -  ka 


(X'.k,  ] 


:*) 


K  I  f2  =  0 


(12a) 


X'Er-X’  -2 
G  ,  p 

E„  +  -E~ 


(i  K*K  +  {  KK*  +  |  k 


2  +  ag < i K*  +  _i_  k 

z  3  A  +  X1  3  A  +  X'  z 


i[M!!Kai^  k*  +  | 

D  X  -3 


+  2  F 


(x',kz] 

X’ 


4h'?zm 


K] 


(x’,kzl 


K)  fl  "  0 


(12b) 


Within  the  framework  of  our  model  (12a)  and  (12b)  arc  exact.  Four  fundamental  energies 
occur  in  Eqs.  (12).  They  are  tabulated  in  Table  3  in  units  of  10  eV  for  InSb  at  experimentally 
interesting  magnetic  field  strengths. 


TABLE  3 

Hydrogenic  Rydberg 
(R) 

Cyclotron  Energy 

(fiw  ) 
c 

eg 

A 

0.6 

10.0  -  100.0 

236.0 

805.0 

It  is  easy  to  estimate  the  order  of  magnitude  of  the  various  terms  appearing  in  (12a,  b)  by 

means  of  the  following  scheme,  valid  for  conduction  band  eigenvalues  and  magnetic  fields  such 

that  hoo  >  R: 
c 
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X'  =  0(1'"  4  fuo  ) 

G  c 


K  or  K  *  =  0[(m*  fiw^)1  2] 


V(r)  =  Q[(Rriojo)1//2] 


(13) 


where  m is  the  conduction  band  mass  in  InSb  (m*  ~  0.0138m^). 

In  view  of  (13),  Table  3  and  the  definition  of  X'  [X'  =  E  —  V(r)),  we  see  that  to  an  excellent 

.  - 1 

approximation  we  can  neglect  commutators  of  K,  Ki,s  and  with  (A  +  X1)  in  (12a,  b). 
Introducing  definitions 

*“E-EG  ■ 


1 


2  1  F  H2  2  n2  2/3(4  +  E) 

Zm*(6)  =  (  3  +  3  A^Te)  E^  honcc  3  |^  =  C/2m*(W  where  C  =  -  -  — 


/  c  \  efl  T?  /  c  \  ^  m  v  ( u ) 

a?  <p)  -  —/  T\ —  t — 5 

c  m  ( to )  c  2  r  2 

n  c> 

o 


Y<&)  =  ficoc(6)/2R(6)  , 


,l/2 


(14) 


and  taking  as  unit  of  length  r^  =  [fi/2m*(&)  w  (&)]  7  and  unit  of  energy  fkj^(  6),  (12a,  b)  become 


H1 11 1  +  II12I2 


H00f  0  +  II 0  .  f 
22  2  21 


(*  +  Eg)  "ft  -  f,f1 

i -hi:)  « 


fofZ  :  fef2 


(15a) 


(15b) 


where,  in  the  gauge  of  section  1,  we  have 


H 


11 


1  +  0.5C 

Cttw 


E  -  V(r  r) 


A  +  E  —  V(r  r)  A  +  E 
c 


)]  (-v2  *?5*!  LJ  *  (*  *  Sr) 


u(r) 


H„  =  H .  .  +  r  C 


1  du(r) 

9 

i  a 

r  dr 

|P  9p" 

‘  2  d(p 

A 

Cfiu 

+  V  +  7--k\ 


_  i  c  ( _ i _ ) 

4  \A  +  E  -  V(rrr)/,  (16a) 


X  du(r) 


22  11  2  A  +  E  —  V(r  r)  E  —  V(  r  r)r  dr 

—  c  c 


i  &-_£ L_ 
top  4 


(16b) 


iC  fico 


II 


12  2(E  —  V(r  r))  r  dr 


t  o  i  o  ,  p  .  .  c;  , 

~z(^  “  p  ^  +  71  +  p  ail 


(16c) 


H21  =  Hi2 


where  u(  r)  =  V(r  r)/fia>  .  Decoupling  (15a,  b)  we  obtain  the  Sehrodingcr-like  equations 

E  \-l 


(nil-H12(H22-fV  H2ll  fl  =  V1 


lH22-HZ1(Hll-Rirll,121  f2  =  V 


(17a) 

(17b) 


6 


Our  desire  to  solve  for  6  and  ft^  from  (17a,  b)  by  the  variational  method  is  thwarted  hy 
the  difficulty  of  evaluating  the  expectation  value  of  terms  containing  the  inverse  operators.  In 
Appendix  A  we  show  how  to  evaluate,  approximately,  the  effect  of  these  terms.  For  the  present 
we  need  only  the  result  that  for  the  eigenvalues  of  the  lowest  lying  M  =  0  or  M  =  1  states  we  can, 
to  very  good  accuracy,  write  simply: 


H,,f.  =  ft  f  (18a) 

111  j  1 

"zz^V2  '  (18b) 


The  reason  for  using  variational  methods  for  estimating  eigenvalues  of  Schrodinger  equations 
is,  of  course,  that  a  relatively  inaccurate  trial  wave  function  can  produce  a  fairly  accurate  en¬ 
ergy.  In  (18a,  b)  the  Hamiltonians  and  are  themselves  functions  of  the  eigenvalue  ft  ; 
moreover,  ft  is  not  itself  an  eigenvalue  but  rather  a  function  of  the  eigenvalue  as  indicated  in 
(15a,  b).  Is  it  true  that  for  (18a,  b)  the  variational  method  produces  a  good  energy  from  an  only 
moderately  good  wave  function?  This  question  is  answered  affirmatively  in  Appendix  B. 

Using  trial  functions  of  the  form  (7a)  and  (7b)  we  have  solved  (18a)  and  (18b)  variationally  by 
computer,  making  sure  that  the  final  variational  energy  deduced  is  equal  to  the  energy  which  was 
inserted  into  H ^  and  H^,  respectively.  Our  results  for  ft  (M  =  0),  the  ionization  energy  asso¬ 
ciated  with  the  spin-up  impurity  ground  state  level,  are  given  in  Table  4.  Values  of  R,  E^,  and 
A  employed  are  taken  from  Table  3.  Energies  are  in  meV. 


TABLE  4 


H(kG) 

12.0 

18.0 

25.0 

50.0 

80.0 

120.0 

$4(  M  =  0) 

1.97 

2.27 

2.54 

3.23 

3.80 

4.40 

Table  4  results  are  obtained  using  a  bottom-of-the-band  Rydberg  of  0.6meV  for  InSb.  Given 
a  different  Rydberg  value,  Hq,  where  0.5  <  Rq  <  0.8,  ft|(M  =  0)  in  meV  can  be  obtained  approxi¬ 
mately  by  multiplying  the  tabulated  value  by  1  +  1.1  (Rq  —  0.6),  whereas  the  analogous  factor  for 
=  1)  is  1  +  (R  -  0.6). 

t  °  T  „  J 

We  have  plotted  the  ionization  energy  difference  ft^(M  =  0)  —  ft^(M  =  1)  =  A^  (H)  in  Fig.  1  both 
for  the  nonparabolic  and  parabolic  bands  with  identical  bottom-of-the-band  mass  and  Rydberg 


The  quantity  A  (H)  is,  in  our  model,  expected  to  be  equal  to  the  energy  difference  between  the 
cyclotron  resonance  and  impurity  cyclotron  resonance  absorption  lines  at  fixed  magnetic  field. 

The  striking  nonparabolic  band  effect  on  A^(H)  shown  in  Fig.  1  can  be  understood  qualitatively 
s  follows.  In  the  InSb  conduction  band  the  electronic  effective  mass  increases  with  electron 
energy.  Thus  the  ionization  energy  of  the  M  =  1  impurity  level  increases  markedly  with  magnetic 
field  relative  to  the  corresponding  nonparabolic  energy.  On  the  other  hand  the  M  =  0  impurity 
level,  lying  always  near  the  bottom  of  the  band,  behaves  much  the  same  for  both  parabolic  and 
nonparabolic  bands.  At  high  fields  the  electron  mass  characterizing  the  M  =  1  level  in  the  non¬ 
parabolic  band  increases  sufficiently  rapidly  that  A^(H)  reaches  a  maximum  and  then  begins  to 
decrease  with  increasing  field.  This  effect  should  be  detectable  experimentally. 

The  spatial  spread  of  the  envelope  function  f^  gives  a  good  idea  of  the  degree  of  localization 
of  the  spin-up  electron  with  respect  to  the  Coulomb  center.  For  the  ground  state  (M  =  0),  we 
find  that  f^  is  somewhat  less  markedly  elongated  in  the  direction  along  the  magnetic  field  than 
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Fig.  1.  Comparison  of  Aj(H)  computed  for  the  nonparabolic  conduction 
band  of  InSb  and  for  a  simple  parabolic  band  with  the  same  bottom-of- 
the-band  parameters,  y^  =  ^/2R  evaluated  at  the  bottom  of  the  band. 


the  wave  function  calculated  by  YKA.  Of  greater  interest  is  the  very  marked  elongation  of  the 

M  =  1  level.  If  wc  let  p  be  the  radial  coordinate  at  which  |f,  |  achieves  its  maximum  value, 

max  1  i  i  2 

|fl  |  ,  at  z  =  0  [that  is,  fj  max  =  fj  (Pmax>  °M,  then  wc  define  the  extent  of  [fj  |  in  the  z- 

direction  by  z  where 
J  o 


|  f  .  (p  ,  z  )/f .  <p 
1  rmax  o  1 


„  0)|  2  =  e'1 


(19) 


Our  calculations  show  that  zq  varies  from  ~2  X  10  ^  cm  at  1  2  kG  to  ~1  x  10  '’at  120kG  in 
InSb.  Since  the  very  purest  samples  of  InSb  presently  available  have  a  mean  interimpurity  sep¬ 
aration  not  significantly  smaller  than  2  X  10"5cm  it  seems  possible  that  impurity  overlap  effects 
may  not  be  entirely  negligible  for  the  M  =  1  state. 

Finally,  we  remark  that,  unlike  the  parabolic  band  case,  the  M  =-l  state  ionization  energy 
is  less  than,  rather  than  equal  to,  the  M  =  1  ionization  energy.  The  energy  differences  are, 
however,  never  greater  than  a  few  tenths  of  a  meV  for  fields  up  to  120kG. 


IV.  POLARON  EFFECTS 

Because  of  the  small  conduction  band  mass  in  InSb,  the  kinetic  energy  of  electrons  bound  to 
donors  in  InSb  exceeds  the  longitudinal  optical  phonon  energy,  fiu  Q,  in  magnetic  fields  of  ~60kG 
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or  more.  One  .night  doubt,  therefore,  that  the  simple  introduction  of  the  static  dielectric  con¬ 
stant  &  into  the  Coulomb  potential  (2)  would  completely  characterize  the  polarization  properties 
of  the  erystal  for  the  whole  range  of  experimentally  interesting  magnetie  fields  (say  0  to  150  kG). 

9 

Platzman  has  shown  on  the  basis  of  the  Frohlich  polaron  model  that  the  proper  treatment  of 
lattiee  polarization  requires  augmenting  the  Hamiltonian  ( 1 )  by  introducing  the  Frohlich  phonon  en- 
ergy  and  electron-phonon  interaction  terms  and  solving  for  the  eigenvalues  of  the  new  Hamiltonian. 
Beeause  we  are  interested  in  relatively  low  lying  states  and  because  the  Frohlich  electron-phonon 
coupling  constant,  a,  is  very  small  (a  —  0.02)  in  InSb,  we  are  justified  in  calculating  corrections 
to  our  impurity  levels  to  lowest  order  in  a  and  in  assuming  a  simple  parabolic  conduction  band 
for  obtaining  these  corrections. 

Thus  for  our  purposes  it  will  suffiee  to  consider  the  Hamiltonian  K  whieh  is  written  below 

using  units  appropriate  to  a  polaron  calculation  (length  in  units  of  r  =  (fi/2mw  )*,  energy  in  units 

of  ftu;  ). 
o 

K  =  K0U)  +  K1  ,  (20a) 

K0(X)  =  7T2  +  U(r)  +  V  b^bk  ,  (20b) 

K1  =  L  ‘V0 - bk  +  e - bk'  '  (20e> 

_L  12 

where  -  (4ra/S)*/k  with  S  =  crystal  volume  and  a  =  y  U/&x  -  i/&0)  e  /r  /fi w 
U(r)  =  XV(ror)/ho?o 

?  i  r  e  1  2 

7r  =  -  V  +  -  A(r  r) I 

l  —  fie  —  o  J 

b ^  is  the  creation  operator  for  a  LO  phonon  of  wave  veetor  k  and  X  is  a  parameter  introduced 
to  facilitate  the  discussion  of  Appendix  C.  For  present  purposes  we  consider  X  =  1 . 

We  will  confine  ourselves  to  obtaining  a  rigorous  upper  bound  to  the  polaron  corrections  to 
the  ionization  energy  of  the  ground  state  of  K  in  lowest  order  in  a.  Except  near  the  polaron 

resonances^  we  would  expect  that  the  shift  in  ionization  energies  of  the  excited  states  would  be 

of  the  same  order  of  magnitude  as  that  of  the  ground  state. 

In  Appendix  C  we  show  that  a  rigorous  upper  bound  to  the  shift  of  the  ground  state  ionization 
energy  due  to  eleetron-LO  phonon  interaction  is  given  to  lowest  order  in  o'  by 

f  fta,o(<0U2|0>-i^)  .  (21) 

o 

where  |  0  >  is  the  ground  state  eigenfunction  of  Hq(\  =  1).  From  (21)  we  ean  eonelude  that  to 
order  a  the  ionization  energy  shift  of  a  polaron  in  the  ground  state  is  overestimated  at  all  fields 

by  replacing  the  band  mass,  m,  by  the  polaron  mass  m*  =  (1  +  a/6)  m  in  (1).  We  believe,  in 
addition,  that  the  replacement  of  m  by  m*  in  (1)  gives  a  quite  accurate  estimate  of  the  polaron 
correction  for  the  range  of  fields  considered  here. 

Numerical  calculation  of  (21)  employing  optimized  wave  functions  of  the  form  (7a)  and  using 
values  of  0.6meV  for  the  hydrogenic  Rydberg,  0.0138  for  the  band  mass  to  vacuum  mass  ratio, 
and  0.02  for  a  shows  that  polaron  effeets  shift  the  ground  state  ionization  energy  certainly  less 
than  0.25%  for  magnetic  fields  from  4  kG  to  120kG. 
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V.  CONCLUSIONS 


Nonparabolic  effects  can  not  be  neglected  in  calculating  the  binding  energies  of  impurity 
states  w  hich  lie  near  or  above  the  n  =  1  Landau  level  in  magnetic  fields  greater  than  ~  1  5  kG  in 
InSb.  The  nonparabolic  effects  result  in  greater  ionization  energies  for  the  higher  lying  states 
than  would  be  expected  from  calculations  for  a  simple  parabolic  band.  These  effects  should  be 
observable  experimentally  in  cyclotron  resonance. 

For  samples  of  highest  attainable  purity,  overlap  effects  between  adjacent  impurity  sites 
should  be  insignificant  for  the  ground  state.  Overlap  may  be  important  for  excited  states,  par¬ 
ticularly  those  which  tend  to  be  strongly  elongated  along  the  magnetic  field  direction  such  as  the 
M  =  ±1  levels. 

Central  cell  corrections,  not  taken  into  account  in  the  present  calculation,  may  be  important 
for  the  ground  state  ionization  energy  and  will  be  expected  to  increase  rapidly  with  magnetic 
field. 

In  the  simplest  possible  picture  we  can  imagine  that  the  central  cell  correction  is  propor¬ 
tional  to  the  probability  of  finding  an  electron  in  a  small  volume  (dimensions  of  the  order  of  the 
lattice  spacing)  around  the  impurity  center.  If  we  assume  that  the  effective  mass  wave  function 

is  contained  in  a  volume  whose  dimensions  are  proportional  to  the  cyclotron  radius,  then  the  in- 

3 /2 

crease  of  the  central  cell  correction  with  magnetic  field  should  be  proportional  to  H  '  . 

For  states  with  M  0  we  would  expect  central  cell  corrections  to  be  relatively  very  small. 

Poloron  effects  on  the  ionization  energies  are  negligible  except  at  fields  near  which  the 
polaron  resonance  occurs. 
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APPENDIX  A 


We  consider  the  energy  shift  introduced  by  the  presence  of  the  perturbations  H12  and  H ^ 
in  (15a)  and  (15b),  respectively.  These  terms  have  been  neglected  in  (18a,  b)  and  the  discussion 
which  follows  in  the  body  of  the  paper. 

First  we  observe  qualitatively  that  from  (16c),  and  are  small  quantities,  being  only 

of  order  (ha;  (  R)*  fiu; c/E.  Further  f^  and  f^,  which  are  coupled  by  the  perturbation,  differ  in  en¬ 
ergy  by  an  amount  of  the  order  of  the  conduction  band  spin  splitting,  which  is  small  compared  to 
E^.  Thus  we  have  a  small  perturbation  coupling  together  levels  which  lie  relatively  near  each 
other  in  energy.  We  do  not  expect,  therefore,  any  significant  shift  in  energy  due  to  and  H^ 
To  make  a  quantitative  calculation  we  focus  for  definiteness  on  (17a),  reproduced  below 


H 


AAiA-  [H.  -  &  )-1H, ,]  {,  =6  f, 

111  1  12  22  f  21  J  1  f  1 


From  ( 1 6b)  we  observe  that  to  a  good  approximation 


H22  H11  +  2  C  A  +  E 


=  H11  +  K 


(Al) 


Inserting  into  (Al)  gives 


’ll  =  H12<H11  +  K-  H 


H 


21 


( A2) 


We  still  cannot  compute  6^  variationally  because  of  the  difficulty  of  evaluating  expectation  values 
of  H'^.  However,  by  exploiting  the  smallness  of  we  can  obtain  accurate  values  for  6^  by 

a  trick.  Since  —  commutes  with  we  can  choose  energy  eigenstates  which  are  simul¬ 
taneously  angular  momentum  eigenstates.  WTe  are  interested  in  finding  the  energies  of  the  lowest 
lying  angular  momentum  eigenstates  with  M  =  — 1,  0,  1.  For  each  value  of  M  we  construct  a  new 
Hamiltonian  whose  eigenvalues  for  given  M  are  equal  to  the  eigenvalues  of  H^  —  for  the 
same  M  to  order  Ih^I^.  However  H^,  unlike  H^  —  H^,  is  so  constructed  that  it  is  easy  to 
find  its  eigenvalues  variationally. 

Such  a  Hamiltonian  is  given  in  (A3) 


hm  =  hm 


+  Hi 


M 


HM  =  HU-R(ctz-1/2) 


Hi 


M 


(PMH12  +  H21PM)  ctx 


(A3) 


where 


vIC-!)  ■  •,-(?}) 


and  P^  is  the  Hermitian  operator  defined  by 

P  eim^=(l-6  M))e^ 

M  m,  M-l 


c(o) 


W7e  easily  verify  that  is  Hermitian,  and  if  Ik  is  an  eigenfunction  of  H^  with  eigenvalue 
then  h . |  t >  and  h.||>  are  eigenfunctions  of  H^  with  eigenvalues  and  +  K.  respectively. 
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_ t  _  2 

Treating  II ^  as  a  perturbation  on  we  find  that  to  order  |H^|  the  perturbed  eigenvalue  of 

Hyj  is  given  in  perturbation  theory  by 

<Htl  ~  1,'11>  •  {A4) 


where  the  expectation  value  is  taken  in  the  eigenstate  of  II ^  of  interest  with  angular  momentum 
eigenvalue  M.  But  (A4)  is  the  expectation  value  we  want  to  find.  Thus  if  we  compute  the  eigen¬ 
values  of  by  an  appropriate  variational  method  we  should  be  able  to  approximate  (A4). 

We  have  list'd  the  following  trial  functions: 

>lcs  =  ag0  |t>  +  ibzg1  |l> 


°  IU  =  1  =  agl  ^  +  iljpz  ^  ?i  l*> 

">  ■  (i)  •  i»  ■ (?)  ■ 


( A5) 


where  gQ,  g^  and  g^  have  the  form  (7a),  (7b)  and  (7b), respectively,  and  both  g^  and  g^  are 
M  =  +1  eigenfunctions  of  L^.  Note  that  neither  ipQ§  nor  are  eigenstates  of  and  that  the 

method  described  here  is  useless  for  approximating  the  eigenfunctions  of  —  H^. 

For  ^Qg  and  we  have  8  and  9  independent  parameters  to  vary  respectively  (b  is  a  vari¬ 

ational  parameter).  In  both  cases  we  find  that  the  effect  of  °n  the  binding  energies  calculated 

from  H.  .  alone  is  less  than  1%. 

1 1 
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APPENDIX  B 


For  definiteness  we  consider  the  error  in  determining  &  by  the  variational  method  for  the 


eigenvalue  problem  (18a) 


'  fl  "  (*  +  EG)t;t  fl 


Vl 


(Bl) 


with  6^  )  given  by  ( 1 6a). 

Assume  that  for  the  exact  solution  of  (Bl),  &exacj.»  1 be  error  in  estimating  the  eigenvalue 
of  H^(&exact>  using  our  variational  trial  function  is  <5 v  where  <5 is  small  enough  to  allow  neglect 
of  quantities  of  order  <5^ .  The  final  result  of  the  variational  calculation  will  be  an  approximate 


value  of  6  6 


approx 


Let 


<5=6  -6 

a  approx  exact 


If  we  can  neglect  terms  of  higher  order  than  <5  then  the  error  introduced  in  the  variational 
calculation  by  inaccuracy  of  6^  in  the  Hamiltonian  )  is 


9HM(6t) 


a  & 


6-6  , 

t  exact 


( B2) 


Thus,  the  total  error  introduced  in  6^  of  (Bl)  is 


6  + 
v 


<»;,> 


■  ( 


so  that 


6  = 
a 


26 

approx 

eg 


<sil> 


exact 


(Rho;  )2 
c 


( B3) 


We  con- 


Note  from  (B2),  (16a)  and  (13)  that  is  a  quantity  of  order 

elude  then  from  (B3)  that  for  (16a)  the  variational  procedure  leads  to  an  estimate  of  6 
whose  accuracy  is  not  worsened  by  the  fact  that  the  eigenvalue  of  H^(6^  )  is  6^  and  not  6  and 
that  itself  depends  upon  6^  .  The  analogous  remark  applies,  of  course,  to 


exact 


\ 
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APPENDIX  C 


To  derive  an  upper  bound  to  the  shift  in  the  ground  state  ionization  energy  due  to  the  presence 
of  from  (20e)  in  the  Hamiltonian  H,  we  treat  as  a  perturbation  on  the  eigenstates  of  HQ(x). 
Denoting  by  |  0,  x}>  and  |n,  X>  the  ground  state  and  n**1  excited  state,  respectively,  of  Kq(X)  and 
bv  Eq^  and  E  the  corresponding  unperturbed  eigenvalues,  we  obtain  the  lowest  order  polaron 
correction  to  the  ground  state  of 


”P01  ~  (2-r)3  U 

'  n  =  0 


v  I  <n,  x|  e1--  |  0,  X>  | 
~EnX~E0x  +  1 


(C'l) 


Following  Platzman  we  iterate  the  identity 


IEn>-E0>-k  > 


F'nX  E0X  +  1  1+k2  (E  -  +  1)  (1  +  k2) 

nX  ux 


(C2) 


to  obtain  from  (CM) 


EpolU)  =  "  T1  +  T2  -  T3W  • 


Mn0(k.X) 


J  1  +  k 


(C3a) 


t,  =  y 

C.  i_i 


I  * 


-E0x-k  ) 


2  2 

(1  +  K  j 


(C3b) 


fivi  . 

dk  — ^ 


M_2A(k.X)(Enx-Eox-k2)2 


2,2 


'**»>  <E„x'Eox*’> 


(t’3c) 


where  M  ^n(k,  X)  =  — v?*  |  <n,  X  |  e*—  —  |  0,  X^  |^.  From  (C 3c)  it  is  important  to  note  that  every 
no  /  *,  \j  K 

(  27t) 

term  in  is  positive. 

fly  completeness  we  have 


T  =  —^3  f  dk  -^L 

1  (  2 7r) J  J  1+1 


Using  the  relation 

V  x  yn\ .-ik-r  I_v  ✓_l_ik*ri,v\  _y„i,.-ik.r  „  ,  .ik-n 

n 


(En-E0)  <0|o"-'-|n>  <n|el--|0>  =  <0 1  (e_1--,  K0]el-^|0>  . 


we  obtain 

T2  =  °  , 

where  we  have  used  <(0,  X  |  k  •  u_  |  0,  x>  =0,  which  follows  from  the  fact  that  I  0,  x^  has  definite 
parity. 
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The  problem  remaining  is  to  approximate  T^(A)  in  an  appropriate 
to  replace  E^  —  EQ^  +  1  in  the  denominator  of  T^(a)  in  (C3c)  by  some 


way.  Our  procedure  is 
funetion  g(k)  thus  obtaining 


x3(x)  =  V 


I  * 


M>-  IEn»  -  Eo>  -  tV 

(1  +  k2)2  g(k) 


=  C(g.  x)  t3(x) 


(C4) 


where  C(g,  A)  is  a  number  which  depends  upon  the  form  of  g  ehosen  as  well  as  A. 

We  now  make  two  very  plausible  assumptions.  First  we  assume  that  if  we  denote 

lim  T-(A)  bv  T^(0+),  then  T-(0+)  =  T-(0)  and  likewise  T^(0+)  =  T^(0).  Second  we  assume  that  if 
A-O-f 

we  denote  lim  E  %  =  E  _  ,  then 
,  nA  n0+ 


En0  +  E00  +  <  Enl  E01 


(C5) 


What  (C5)  says  is  that  turning  on  the  Coulomb  potential  lowfers  the  ground  state  energy  more  than 
it  lowers  the  energy  of  any  other  level.  This  is  obviously  true  w'hen  there  is  no  magnetie  field. 

It  is  plausible  generally  because  the  ground  state  is  the  only  nodeless  state  and  therefore  can 
crowrd  itself  in  toward  the  Coulomb  center  with  minimum  cost  in  eurvature  (kinetic  energy)  com¬ 
pared  to  any  other  state. 

From  (C5),  (C4)  and  (C3e)  wre  eonclude  that 

C(gfl)  >  C(g,  0+)  .  (C6) 

The  polaron  correction  to  the  ground  state  ionization  energy  is  by  definition 

Epol(°+)“Epol(1)  =  T3(l)-T3(0)  •  (C7) 

The  quantity  T^fl)  —  T^(0)  will  certainly  be  an  upper  bound  to  T^fl)  —  T^(0)  if 


C(g,  0)>1  (C8) 

by  virtue  of  (C4)  and  (C6). 

We  ean  ensure  the  validity  of  (CS)  by  setting  g(k)  =  1  as  is  obvious  from  (C3c)  and  (C4).  In 
that  case  we  evaluate  (C4)  using 


^(En-E0)2  <0 1  e_i-  -  |  n>  <n\e^\0> 

n 

=  <0|  f[e_i— Kq],  K0}  e1— ”-|  0> 

=  <0  |  ( 2k  •  I  +  k2)2  I  0> 

and 


[e"^E  M01  =  (2k  .  w  +  k2)  e_i—  - 

to  obtain 

T3(l)  -  T3(0)  =  |  a  (<1,  0  I  7r2  I  0,  1>  — 
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Computer  evaluation  of  (Cl)  at  \  =  0  as  a  function  of  magnetic  field  and  asymptotic  forms  for 

It  2 

small  and  large'  field  show  that  (C8)  is  also  satisfied  by  setting  g(k)  =  1  +  k  .  With  this  choice 

of  g  \vc  find 

T3(l)  -  T3(0)  <  T3(l)  -  T3(0)  =  2.  [<1,0  I  r2  |  0,  1>  -  .  (C9) 

\  u  o  ' 

Finally,  we  remark  that  setting  g(k)  =  1  +  k^  is  equivalent  to  applying  (C2)  to  (C3c)  and 
truncating  the  resulting  series  for  T^(A)  after  the  first  term.  This  procedure  generates  the 
polaron  effective  mass  approximation  for  the  problem  of  a  polaron  interacting  with  a  central  po¬ 
tential  and  applied  magnetic  field.  For  weak  magnetic  and  Coulomb  fields  the  polaron  effective 
mass  approximation  and  hence  our  bound  (C9)  should  be  excellent. 
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